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COMPUTER-AIDED ANALYSIS OF LANDSAT DATA 


FOR SURVEYING TEXAS COASTAL ZONE ENVIRONMENTS* 

S. J. Kristof and R. A. Welsmlller 

ABSTRACT 

A study was conducted to detennlne the feasibility of using machine- 
aided processing of Landsat data to Inventory environmental units within 
the Texas coastal zone. The analysis was conducted on geometrically 
corrected and spatially registered Landsat data collected on November 27, 
1972 and February 25, 1975 over the Matagorda Bay area of the Texas coastal 
estuarine system. A clustering algorithm (nonsupervised processor) was 
used to divide the data into groups of sample points of similar spectral 
characteristics. Correlation of spectral classes with reference data on 
a point-to-point basis showed the coastal features exhibit unique spectral 
variations. Statistics developed on these groupings were input to a maxi- 
mum likelihood algorithm and the test sites classified. The following 
terrestrial and aquatic environments were discriminated: alternating 

beach ridges, swales, sand dunes, beach birms, deflation surfaces, land- 
water Interface, urban, spoil areas, fresh and salt water marshes, grass 
and woodland, recently burned or grazed areas, submerged vegetation and 
waterways . 


*This work was supported by the National Aeronautics and Space Adminis- 
tration under contract NAS9-14016 and NAS9-14970. Journal Paper Number 
6890 Purdue University, Agricultural Experiment Station. 
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Visual observation of results obtained from both the November 1972 and 
February 1975 data Indicate that no major differences existed in the land 
resource maps* The results did show that analysis of Landsat data with computer- 
aided techniques is a viable technique for surveying coastal features. 
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INTRODUCTION 

Land resources Inventories of coastal environments are vital to the 
prudent management of these resources by national, state and local govern- 
ments. Remote sensing techniques provide a means of obtaining this neces- 
sary information (2). Vytautas, et al. (6) found multispectral analysis 
of aerial imagery to be a valuable tool in the mapping of the Delaware 
wetlands. Using machine— aided processing of Landsat data, land resources 
inventories of large areas can be produced with a minimum of time and 
expenditure. Tamocai and Kristof (5) used this technique to survey the 
l^ckenzie River delta area in Canada. This paper examines the feasibility 
of using machine-aided processing of Landsat data to inventory pertinent 
environmental units within the Matagorda Bay region of the Texas coastal 
zone . 


DATA 


Landsat multispectral scanner (MSS) data collected on November 27, 

1972 and February 25, 1975 were used for this study. These data, covering 
the Matagorda area (Figure 1) , were geometrically corrected (1) , overlaid 
one to the other and spatially registered to ground control points selected 
from U.S. Geological Survey (USGS) 7% minute topographic quadrangle maps 
(Figure 2) . The resulting multichannel, multidate data set when printed 
on a computer line printer in pictorial form had a scale of 1cm = 2400cm 
and could be physically overlaid on USGS topographic maps of the same 
scale. Thus, individual pixels could be located for accurate ground 




observation. 
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STUDY AREA 

Two study areas in the Matagorda Bay region of the Texas coastal 
zone were selected for this investigation. These areas are represented 
by the Pass Cavallo and Port O’Connor DSGS 7% minute quadrangles. The 
Pass Cavallo area consists of Matagorda Island and several other small 
barrier Islands. The dominant geomorphlc features of the area are the 
alternating ridges and swales paralleling the beach of the gulf side. 

Since the ridges retard runoff, the swales tend to be moist and period- 
ically Impound water. The Port O'Connor area, situated between the 
coastal wetlands and the Inland woodlands, consists chiefly of sand 
and clay material deposited from overbanking streams and sediments 
during the Pleistocene history of the area. It is traversed by elongated 
sand belts with very slight topographic relief. In general, the entire 
area appears as a broad, coarse textured terrain with fresh-water marshes 
and oak motte environments. 

Both areas. Pass Cavallo and Port O’Connor, represent an interactive 
ecological region and are frequently affected by wind tides, tidal creeks, 
shoreline erosion, and active surface faulting. The native vegetation 
stands represent a collection of vegetation of similar species dominated 
by a single characteristic species of a small number of co-dominated 
species. The lower areas are covered with marshes and swamps. Salt-water 
marshes are characterized by pure stands of Spartlna alteniflora (cord- 
grass) , Salicomia perenls (glasswort) and Suaeda spp. (seepweed) . Closed 
brackish-water marshes occur in the central inland area of Matagorda 
Island and other barrier Islands and on the shoreface of the mainland 
sides of the bays. Periodic salt-water inundation and fresh-water runoff 
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Figure 2 


Portions of the Port O’Connor and Pass Cavallo, U.S. Geological Survey 7^5 minute topographic 
quadrangle maps. 
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from adjacent higher lands constrain biologic growth In these areas. 
Vegetation primarily consists of Spartlna patens (marsh hay cordgrass) , 
Spartlna cynosuroldes (big cordgrass), Sclrpus spp. (bullrush) , Typha 
latlfolla (cattail), and Juncus spp. (rushes). At higher elevations 
fresh-water marshes are developed and maintained by rainfall or a perma- 
nently high water table. Fresh-water marshes are characterized by Juncus 
spp. (rushes), Sclrpus spp. (bullrush), Typha latlfolla (cattail), and 
Spartlna pectlnata (sloughgrass) . 

Topographically low frequently Inundated areas may or may not be 
occupied by vegetation. Aquatic vegetation Is followed by zones with 
Sallcomla spp. (glasswort) and Dlstlchlls splcata (spike grass). The 
drier parts of the marsh are characterized by Spartlna patens (marsh hay 
cordgrass) and Spartlna spartlnae (coastal sacahulsta) . Swamps are com- 
posed entirely of fresh water and are characterized by heavy growth of 
Quercus nigra (water oak), Nyssa blflora (gum), Ulmus spp. (elm), Vltls 
spp. (grapevine), and Mbrun spp. (mulberry). 

The flat uplands which occur landward from the bays and lagoons and 
extend from sea level to an elevation of approximately 100 feet above MSL 
are not significantly affected by water Intrusion. Their main herbaceous 
and woody vegetations are Andropogon spp. (bluestem), Sorghastrum spp. 
(Indlangrass) , Paspalum spp., Prosopls spp. (mesqulte). Sorghum halepense 
(Johnson grass), Celtls spp. (hackberry) , Acacia fameslana (hulsache) , 
chaparral, cactus, and Quercus Virginia (live oak) (3,4,8). 
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METHOD 

Single date, unsupervlsed classifications of the Pass Cavallo and 
Port O’Connor areas using the registered data from November 1972 and 
February 1975 were produced for comparison to each other. False color 
Images (Figures 3 through 10 , presented In black and white) produced 
from a digital display unit were used to obtain an overview of the two 
study areas and to obtain some detailed Information about the surface 
features. The training areas were selected from the Images to contain 
typical examples of each cover type of Interest. 

For each date, cluster analyses were conducted. Three training 
sites for both the land and water areas for each study area were analyzed. 
The training areas for the Port O'Connor area are shown In Figure 11. 

For each study area twenty-one spectrally separable classes resulted from 
the clustering process. Fourteen and fifteen land and seven and six water 
classes were Identified for the Pass Cavallo and Port O'Connor areas, 
respectively. The cluster output Included a cluster map showing the loca- 
tion of the spectral classes and a punched output of field description 
cards for each of the cluster classes. The 4 field description cards 
were input Into a statistics processor to obtain the mean spectral response 
of each cluster class In all four Landsat channels and their covariance 
matrices. By summing the mean relative spectral response values of all 
four bands, the magnitude of the responses for each class can be obtained. 
The mean spectral response values were also used to calculate a ratio A >■ yj 
By observing the magnitude M and the ratio A, the cluster classes can be 
tentatively Identified. 

*V Is the relative Intensity of the visible wavelengths (0.5-0.6ym) + 
(0.6-0. 7ym) and IR Is the relative Intensity of the reflective Infrared 
wavelengths (0.7-0.8ym) + (0.8-1. 1pm). 
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The covariance matrix and mean vector statistics of these classes 
were Input Into a maximum likelihood algorithm which classified the areas 
under Investigation. The classification results were output In alpha- 
numeric form using a threshold value of 0.5. Thresholding delineated 
those points not adequately represented hy the training statistics. New 
training fields were selected from the thresholded areas and their cova- 
riance matrix, mean vector statistics and ratios calculated. The statistics 
of these new classes were merged with the previously calculated statistics 
and the study areas reclassified. A flow chart representing this analysis 
technique Is shown In Figure 12. 

Examination of aerial photographs , field observations , USGS topo- 
graphic maps and Bureau of Economic Geology (BEG) land use and bathymetry 
maps and "Spectral Environmental Classifications" furnished by Lockheed 
Electronics Company (LEG) /Johnson Space Center, confirmed the Identification 
of the classes as proposed by the magnitude and ratlolng Information. These 
same materials were later used for evaluation of the classification results. 

RESULTS 

Computer-derived classifications based upon the separation of surface 
features with different spectral responses are presented as land resources 
maps for the Pass Cavallo and Port O'Connor areas (Figures 13 through 18). 
The Pass Cavallo area Is dominated by Matagorda Island and the surrounding 
coastal waters. The Island Is characterized by sand, brush, marshlands, 
and an abandoned Air Force Installation. Elevation of the Island ranges 
between 0-2 meters. The dominant geomorphlc features are the alternating 
ridges and swales, sand dunes and tidal flats on the lagoon side of the 
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Figure 11. Six training sites selected for 
the Port O'Connor study area. 


Island. Tabular classification results (Table 1) Indicate that the Pass 
Cavallo region Is composed of 80% water, 6.4% swales and ridges, 2.1% 
burned and flooded, 2.3% sand and sand dunes, 1.3% housing and runways, 

6.9% salt-water marshes, and 0.8% non-vegetated areas. 

The Port O'Connor region Is characterized by broad areas of relatively 
flat coastal plains occurring Inland from the coastal marshes, bounded by 
stream dlsected wooded areas on the north. Extensive marsh-covered areas 
at 1.5 m or less above sea level stretch along Matagorda Bay proper. 
Tabular classification results (Table 2) show that the Port O'Connor area 
Is composed of 50.3% water, 9.3% fresh-water marshes, 6.3% salt-water 
marshes, 9.9 grasses, 4.3% swamp, 4.1% burned land, 5.7% shrub and trees. 
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Figure 3. 


Image of Pass 
Cavallo; produced on 
electronic printer/ 
plotter from Landsat 
MSS band 0.5-0.6ym 
data - February 1975. 



Figure 4. Image of Pass 

Cavallo; produced on 
electronic printer/ 
plotter from Landsat 
MSS band 0.6-0. 7ym 
data - February 1975. 
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Figure 5. 


Figure 6. 


Image of Pass 
Cavallo; produced on 
electronic printer/ 
plotter from Landsat 
MSS band 0.7-0. 8pm 
data - February 1975 




Image of Pass 
Cavallo; produced on 
electronic printer/ 
plotter from Landsat 
MSS band 0.8-1. 1pm 
data - February 1975 
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Figure 7. 


1 

i 


Image of Port 
0 ’ Connor ; produced 
on electronic printer/ 
plotter from Landsat 
MSS band 0.5-0.6ym 
data “ February 1975. 



Figure 8. 


Image of Port 
0^ Connor; produced on 
electronic printer/ 
plotter from Landsat 
MSS band 0.6-0. 7ym 
data - February 1975. 
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Figure 9. 


Image of Port 
O'Connor; produced 
on electronic printer/ 
plotter from Landsat 
MSS band 0.7-0.8ym 
data - February 1975. 



Figure 10. 


Image of Port 
O'Connor; produced 
on electronic printer/ 
plotter from Landsat 
MSS band 0.8-l.lvm 
data - February 1975. 





- 15 - 



Figure 12. Analysis flowchart for land use determination. 




Figure 13. Level I Computer Land Resources Map of 
Pass Cavallo (Landsat data collected 
February 1975) . 
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Figure 14. 


Level III Computer Land Resources Map of Pass Cavallo ' 
(Landsat data collected' February 1975) . 
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Figure 15. Level III Computer Land Resources Map of Pass Cavallo 
(Landsat data collected November 1972) . 
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Figure 16. Level I Land Resources Map of Port O’Connor 
(Landsat data collected February 1975) . 
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Figure 17. 
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(Landsat data collected February 1975). 

Vegetated - dry - urban 

Vegetated - dry - mixed - prairie/ cultivated grass environment 

Vegetated - wet - herbaceous - emergent - grazed 

Vegetated - dry - woody/herbaceous 

Inundated - swamps 

Vegetated - dry - prairie grasses 

Vegetated - fresh water marsh 

Vegetated - wet - shrub/trees dominant environment 
Vegetated - fresh water marsh 
Bare land - wet - burned 
Bare, very wet - burned 
Beach sand - dry 
Sand - sand bars 
Salt water marsh 
Bare soil - wet 

’If: 


■ 

' * ; 






21 - 



Figure 18. Level III Land Resources Map of Port O'Connor 
(Landsat data collected November 1972) . 
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Table 1. Tabular output for the Pass 


Group 

Acres 

Beach sand 

464.2 

Sand dunes 

171.6 

Runways 

371.8 

Swales 

1724.8 

Ridges 

925.1 

Non-vegetated 

346.5 

Burned 

576.4 

Flooded 

306.9 

Muddy - Flat 

2180.2 

Submerged (Marshes) 

680.9 

Sand - Shell 

324.5 

Housing 

191.4 

Water 

32976.9 


Cavallo area. 

(February 1975) 

Hectares 

Percentage 

187.9 

1.12 

.69.5 

0.43 

150.5 

0.90 

698.3- 

4.18 

374.5 

2.25 

140.3 

0.84 

233.4 

1.41 

124.3 

0.74 

882.7 

5.28 

275.7 

1.65 

131.4 

0.78 

77.5 

0.46 

13351.0 

79.96 


Total 


41241.2 


16696.8 


100.00 
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Table 2. Tabular output 

for the Port 

O'Connor area. 

(February 1975) 

Group 

Acres 

Hectares 

Percentage 

Urban area 

325.6 

131.8 

0.74 

Prairie grasses 

4314.1 

1746.6 

9.87 

Woody /herbaceous 

1977.8 

800.7 

4.53 

Emergent flora 

1339.8 

542.4 

3.05 

Swamp 

1879.9 

761.1 

4.31 

Fresh-water marsh 

4042.4 

1636.6 

9.25 

Salt-water marsh 

2725.7 

1103.5 

6.25 

Shrub /trees 

526.9 

213.3 

1.21 

Burned vegetation 

1784.2 

722.3 

4.08 

Beach sand 

789.8 

319.8 

1.80 

Sand 

1486.1 

601.6 

3.48 

Bare soil - wet 

512.6 

207.5 

1.16 

Water 

21937.9 

8881.7 

50.27 

Total 

43642.6 

17669.1 

100.00 
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0.7% urban, 6.4% sand and bare soil, and 3.1% emergent flora. 

Classification performance was determined by comparing each classi- 
fication results map with several sets of reference data at a common 
scale, especially the "Spectral Environmental Classification" overlays. 

A hierarchy of land resources features (Figures 19,20) established for 
this investigation, based upon statistical data obtained from LARSYS 
processors and reference data, shows that the 21 classes for each date 
inventoried can be divided Into three distinguishable levels. The first 
level consists of three major classes comprising 22.3% and 3.8% land, 

27.5% and 16.3% wetland, and 50.3% and 80.0% water for the Port O'Connor 
and Pass Cavallo areas, respectively. The lan^ wetland and water group- 
ings can be further subdivided into Level II and Level III classes as 
shown In Figures 19 and 20. 

The land category was further divided into vegetated regions, 
where the vegetation Is not significantly affected by water such as 
native and cultivated grasses; vegetative cover in housing areas; and 
bare, dry areas such as sand dunes, sand beaches and man-made objects. 

The wetland category, where the vegetation Is significantly affected by 
being intermittently Inundated with water, is divided into swamps, marshes, 
swales, bare land covered with shallow water, wet non-vegetated areas, 
and submerged vegetation. 

Water categories were divided into spectral classes based upon their 
chemical and physical properties. Figures 21 and 22 show the spectrally 
separable classes of water for these two areas. These results are In 
agreement with the bathymetry maps prepared by the Bureau of Economic 
Geology, State of Texas. 
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(Landsat data collected February 1975) 
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Figure 21. Spectral Classification of Water for Pass Cavallo Quadrangle 
(Landsat data collected February 1975) ♦ 
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22. Spectral Classification of Water for Port O’Connor 
Quadrangle (Landsat data collected February 1975). 
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At Level III, the rangeland category was separated into woody (oak 
mottes) , herbaceous (native and cultivated prairie grasses) and a mixture 
of woody /herbaceous . The marsh category is also divided into salt-water 
marshes and fresh water and both fresh and salt-water areas where vegeta- 
tion is present. The non-vegetated category is separated at Level III 
into burned, flooded and wet bare land. 

These spectral classes result from differences in the type of yege- 
tative cover, the density of the vegetation, the wetness in the terrestrial 
ecosystems, geomorphic features, and the turbidity, amount of sediment and 
water depth in the aquatic ecosystems. 

Results from the two sets of data, November 1972 and February 1975, 
show little differences in the type and detail of surface features that 
could be delineated. However, results from the November 1972 data show 
the area of beaches and sand dunes to be slightly larger, and burned areas 
were present in the February 1975 data which did not occur in the November 
1972 data. 


CONCLUSIONS 

The results of this investigation indicate that digital analysis of 
Landsat data can be a valuable tool for inventorying environmental resources 
in coastal zone areas. The unique spectral responses of surface features 
within the Texas coastal zone allowed for easy separation of land, water, 
and wetland categories. Classes in vegetated areas were further separated 
based upon the density of grotjnd cover. However, in marsh areas the reflec- 
tance from the background of water and wet soil played a major role in 
class separation. The spatial resolution of the Landsat scanner was 
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adequate to recognize gross patterns of the vegetation, but was not ade- 
quate to resolve individual vegetation species. The spectral separability 
of water classes was a ftmction of turbidity, depth and other physical and 
chemical properties. A visual comparison of results from data collected 
on November 27, 1972 and February 25, 1975 indicated that there were no 
major differences in classification performance between the two data sets. 
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